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ABSTRACT 


AERODYNAMIC ANALYSIS 


A preliminary design tool for advanced pro- 
pellers has been developed combining a fast 
vortex lattice aerodynamic analysis, a fast sub- 
sonic point source noise analysis, and an opti- 
mization scheme using a conjugate directions 
method. Twist, chord and sweep distributions 
are optimized to simultaneously improve both the 
aerodynamic performance and the noise observed 
at a fixed relative position. 

The optimal noise/performance trade-offs 
for straight and advanced concept blades are 
presented. The techniques used include increas- 
ing the blade number, blade sweep, reducing the 
rotational speed, shifting the spanwise loading 
and diameter changes. Some methods yield im- 
provement in both efficiency and noise. 


INTRODUCTION 

PREVIOUS WORK Q,2)* HAS IDENTIFIED THE per- 
formance gains and noise reductions available 
through the use of advanced concept propeller 
blades. These works optimized the performance 
only; noise was considered as a separate issue. 
To properly investigate advanced propeller con- 
cepts as well as conventional geometries, a 
design tool should be able to optimize both per- 
formance and noise at the same time. This paper 
describes the development of such a design tool 
and the noise/performance trade-offs available 
from optimal geometries. 

The design program consists of three parts: 
(1) a vortex lattice method for the aerodynamic 
performance prediction; (2) a subsonic compact 
source model in the time domain for the noise 
prediction; and (3) an optimization scheme using 
a conjugate directions method. 


^Numbers in parentheses designate references at 
end of paper. 


The vortex lattice method has been used by 
various researchers ( 2 to 1) in the study of 
propeller performance. The basic calculation 
geometry is shown in Fig. 1. There is a single 
curved lifting line divided spanwise into a set 
of horseshoe vortices. This lifting line is 
along the quarter-chord line of the blade. At 
control points on the three-quarter chord line, 
a boundary condition of no flow through the mean 
camber surface is imposed. The trailing vor- 
tices representing the wake are prescribed to 
lie on the surface swept out by the quarter- 
chord line with pitch defined by the advance 
ratio. This method is a steady symmetric anal- 
ysis, so the radial circulation distribution is 
the same for all blades. 

The equations to be solved are based on the 
boundary condition that there is no flow through 
the chord plane surface at the control points. 
The boundary condition for the j C ^ vortex 
lattice element is therefore: 
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C ^ 0 zero lift drag coefficient 

design 

C* design lift coefficient 

c chord length 

K drag polar constant. 

As a verification, Figs. 2 to 5 compare the 

predicted results with experimental test data. 
The data in Figs. 2 and 3 are for a representa- 
tive general aviation propeller tested at the 
NASA Lewis Research Center (_ 8 ). The data in 
Figs. 4 and 5, taken from Ref. (^) , are for a 
quarter scale general aviation propeller modeled 
after the blade on a Cessna 172: the McCauley 

1C160. In both cases, the predicted performance 
agrees well, but the predicted variation of 
power coefficient ( Cp ) with advance ratio (J) 
has more curvature tnan the experimental 
results. The fixed wake model is partially to 
blame: the induced velocities would increase 

the wake helix angle, and hence would increase 
the induced drag at the higher loadings. This 
factor would raise the high loading prediction 
of Cp and straighten the curve. 

In these figures, the following symbols are used 

B blade number 

£ 3/4 blade pitch at the 3/4 radius 

J V/nD = advance ratio 

V axial inflow velocity (flight velocity) 

X V/j> 2 k = J/tt = alternate advance ratio 

C-j- T/(p 0 n^D^) = thrust coefficient 

Cp P/(p c n^D^) - power coefficient. 

NOISE ANALYSIS 

The noise analysis method is based on the 
work of Lowson ( 10) . This analysis assumes 
that there is no mean flow (static operation) 
and that the continuous distribution of noise 
sources on the blade can be represented by point 
sources. These assumptions are met when consid- 
ering low speed operation, as in take-off and 
climb, such that the helical tip speeds are less 
than about Mach 0.73. 

Lowson did not consider the noise produced 
by a radial force in his analysis. So, an addi- 
tional term for the radial force noise has been 
derived for use in this study. The equation 
which must be solved for the normalized pressure 
perturbation produced by the propeller is: 
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For the noise analysis, the x-axis is the for- 
ward axis of rotation (m). This equation is 
evaluated at the normalized retarded time 9, 
for the observer time determined by: 



where 

A 0 s sweep introduced delay. 

The propeller is modeled as a distribution 
of point noise sources with strengths propor- 
tional to the local loading and thickness. Each 
of the vortex lattice elements is modeled as two 
noise elements in the spanwise direction (or 
along the hub to hub quarter-chord line for the 
bi-blade) to improve the assumption of compact- 
ness. The layout of noise sources is shown in 
Fig. 6. 

The force coefficients, ( AC^, ACq ,ACr) , on 
the noise elements are half of the coefficient 
on the parent aerodynamic element, and the 
thickness source term is written in the form of 
a force coefficient (AKq) for calculation ease. 
The strength and positions of the source - sink 
pair used to represent thickness are chosen so 
that the Rankine body produced has the same 
chord as the local blade section, and the same 
frontal area as the noise element. 

To increase the execution speed, the pres- 
sure perturbation (Eq. (6)) is calculated at 
equally spaced points in retarded time, 6, along 
with the associated observer times, t. These 
pressures are then interpolated to an equally 
spaced set in observer time. The advantage of 
this method, due to Munro ( JjO , is that the time 
consuming solution for the retarded time associ- 
ated with a given observer time is eliminated. 
Also, the distribution of points in observer 
time is densest in the region of the pressure 
peaks which results in an accurate interpolation 
of the peaks. The multiple blade waveform is 
created from this single blade waveform by 
superposition with the proper phase delays. 

Once the perturbation pressure time history 
has been produced, a Fast Fourier Transform is 
used to obtain the harmonic magnitudes. The 
overall sound pressure level (OASPL) noise meas- 
ure is then calculated from these harmonics both 
in dB and using the A-weighted scale for dBA. 

Figure 7 compares the predicted and meas- 
ured near field noise, one tip radius away in 
the plane of rotation for the quarter scale 
propeller test data in Refs. (9) and (12) to 
(J^). There is a phase shift between the pre- 
diction and the experimental data, but the wave- 
form shape and magnitude are in excellent agree- 
ment. Additional data is available both ahead 
of and behind the plane of rotation. These data 
are compared with predictions in Ref. (15). The 
comparisons show the same phase shift along with 
a small magnitude difference. This magnitude 
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difference results in a few tenths of a dB dif- 
ference in the OASPL. 

OPTIMIZATION METHOD 


r(i) blade element circulations 

a the relative importance of efficiency 

versus noise 


An optimization method automates the proc- 
ess of choosing values of the independent vari- 
ables to create a "best" result. In this work, 
the independent variables are the twist, chord, 
and sweep distributions, while the measure of 
the result is called the Index of Performance. 
The Index of Performance (IP) contains both the 
aerodynamic performance and noise of the pro- 
peller and is formulated so that an optimum 
geometry corresponds to a minimum of the func- 
tion. Different optimal geometries are obtained 
by changing the relative importance of noise 
versus efficiency in the IP. 

The conjugate directions method used for 
the optimization is that of Powell (16) as modi- 
fied by Zangwill ( 17) . After searching along 
each of the current set of search directions, 
the technique searches along the direction re- 
sulting from the end points of the intermediate 
steps. This composite direction is the candi- 
date conjugate direction. If the search proves 
fruitful, the conjugate direction replaces one 
of the previous search directions. 

The line searches are done using a quad- 
ratic fit to three function evaluations. If the 
estimated minimum lies outside the fit interval, 
the interval is shifted towards the minimum 
until the minimum is bracketed. Then, if the 
prediction for the value of the objective func- 
tion does not agree with the computed value, the 
fit interval width may be reduced. The minimum 
interval width allowable is based on the rela- 
tive magnitude of the random portion of the 
objective function. This method is rapid and 
provides some immunity to numerical noise in the 
objective function. 

The Index of Performance (IP) incorporates 
the constraints in the manner of the SUMT form- 
ulation developed by Fiacco and McCormick (18) : 
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C-j propeller thrust coefficient 

Cp propeller power coefficient 

(-target parameter to obtain the desired C p 

N OASPL noise measure in (dB) or (dBA) 

^limit limit value of total sweep 
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tangent to the bound segment between 
nodes j and j + 1. 


Because of the numerical noise present in 
the calculation, the weightings in (8) are not 
decreased to zero as would be customary. In- 
stead, the weightings are held constant and the 
offset target values are adjusted. This pre- 
vents the amplification of the random noise in 
the IP. 

To further speed the optimization, the 
scheme is closely coupled with the aerodynamic 
analysis routine. This is done by recalculating 
only those portions of the model affected by 
geometry changes. For example, if an element is 
retwisted, only the induced velocities at the 
control point change. The induced velocities at 
the bound segment and all induced velocities on 
the other elements remain the same and are re- 
used to calculate the performance. This is most 
beneficial initially when the search directions 
correspond to single geometric quantities. Com- 
pared with a complete solution for a design 
state, this technique is 2 to 24 times faster 
for evaluations of a 12 element blade model. 

An initial test case for the optimization 
method was to maximize straight blade ideal 


4 



efficiency as a function of the twist distribu- 
tion. Goldstein's circulation distribution 
( 19) , which yields the maximum ideal efficiency 
under these assumptions, was obtained, verifying 
the method. 

For a more realistic prediction, a viscous 
drag estimate is included. This estimate 
assumes a laminar flow section operating in the 
laminar bucket with the same dependence of 
on section thickness as for a Clark-Y airfoil 
(to be conservative). The C& for stall was 
chosen to be 1.2, and stall was modeled by a 
sharp increase in the section drag coefficient. 
This drag model is justified because given the 
design lift coefficient, the corresponding air- 
foil shape can be found. By including viscous 
effects, the chord distribution can be a design 
variable in addition to the twist distribution. 

To investigate the trade-offs involved when 
noise is considered, various techniques for noise 
reduction are applied to a straight, a propleted 
and a bi-bladed propeller. These techniques 
include: reducing the diameter to reduce tip 
speed, shifting the loading inboard, and reduc- 
ing the rotational rate to reduce tip speed. 

For these trade-offs, the design variables are 
the twist, chord and sweep distributions. 

Noise reductions are achieved through two 
mechanisms: phase cancellation and noise source 
strength reduction. Phase cancellation is 
accomplished either by adding blades or through 
sweep. The source strength dependence on pro- 
peller parameters depends on the type of noise 
source. The pressure perturbations due to 
thickness terms are proportional to relative 
Mach number squared, and the perturbations due 
to force terms are proportional to Mach number 
cubed. Loading noise is linearly proportional 
to the blade loadings though. Hence, reducing 
the relative Mach number through a reduction of 
the rotational speed or a reduction of the pro- 
peller diameter are expected to have a large 
effect on the noise. 

The resulting optimized trade-offs are shown 
in Figs. 8 to 10 for the three blade types. The 
noise is calculated at a fixed position, 333 tip 
radii away, and 30° behind the plane of rota- 
tion. This angle is near the maximum of the 
directivity pattern, and at a distance appropri- 
ate for sideline measurements in take-off. The 
trade-off curves were obtained by varying the 
parameters in the IP, Eq. (8). The weighting 
value, a, was shifted from consideration of 
only aerodynamic performance towards considera- 
tion of noise only. For each value of the 
weighting parameter, the propeller geometry was 
re-optimized. These optimal geometries produce 
the curves shown. 

Figure 8 shows the trade-off results for 
the straight blade. A simple and very effective 
technique for performance improvement is to 
increase the blade number. The blade number as 
shown starts at three and increases in incre- 
ments of one. When either the propeller solid- 
ity is held constant, or the blade planform is 
held constant, the noise is decreased while the 
efficiency goes up. Additional blades reduce 


the noise by introducing phase cancellation for 
certain harmonics while the strengths of the 
remaining harmonics are constant. The effi- 
ciency is increased because of reduced tip and 
swirl losses. The constant solidity cases have 
higher efficiencies than the constant planform 
cases because the viscous drag is lower. 

Sweeping the blade is the next most effec- 
tive method shown. This trade-off starts with 
no sweep, and ends up with the tip swept back 
about 45° from the initial geometry. The points 
along the curve represent equal increments in 
the weightings, but unequal changes in sweep. 
Sweep reduces the noise by introducing a phase 
shift and subsequent phase cancellation between 
noise from different sections of the blade. If 
the swept blade is concurrently retwisted to 
maintain the radial load distribution, the effi- 
ciency will remain nearly the same. A drawback 
to sweep is that for a blade with a low rota- 
tional Mach number, the amount of sweep needed 
for a useful reduction in noise is large. 

Reducing the rotational speed can reduce 
the noise with a small efficiency penalty. The 
changes shown correspond to 5 percent changes 
in rotational speed. The higher advance ratio 
results in higher blade loadings, hence greater 
tip and swirl losses, but there is a large 
decrease in the noise produced. 

Shifting the loading inboard can reduce 
the noise, but with a sizeable efficiency pen- 
alty. The efficiency drops because the effec- 
tive diameter is reduced, increasing the disk 
loading. The inboard shift of the loading 
reduces the noise by reducing the strength of 
the high Mach number sources, and increasing 
the strength of the low Mach number sources. 
Although hard to quantify, the loadings shown 
are shifted so that the peak moves from the 
75 percent radius inboard to about the 
60 percent radius. 

A very ineffective method for reducing the 
noise is shown to be reducing the diameter. 
Changes here correspond to 5 percent reductions 
in diameter. The noise is reduced by the 
removal of the high Mach number tip sources. 

The remaining sources are increased in strength 
though, which offsets some of the reductions. 

The efficiency drops sharply because of the 
increased disk loading. This performance pen- 
alty is large and makes the method undesirable. 

Figure 9 shows the trade-off results for a 
propleted blade. The character of some of the 
trade-off curves has changed from the straight 
blade results. Most notable is the blade number 
trade-off which no longer shows an increase in 
efficiency with the decrease in noise. This 
curve is for a constant planform; for structural 
reasons, the tip thickness and chord are fixed. 
With increasing blade number, the constant blade 
planform propeller has increased viscous drag. 
Since the tip device is most beneficial at 
higher blade loadings, increasing the blade 
number will also reduce the efficiency gain. 

This reduced efficiency gain, along with the 
viscous drag increase, combine to reduce the 
efficiency of a propleted propeller. Further, 
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while Che total loading remains approximately 
constant, the additional tip thickness noise 
offsets the noise reductions from increased 
blade number. 

Loading shift and diameter reduction show 
unfavorable trends similar to those for the 
straight blade. Loading shift penalizes a pro- 
pleted propeller more by unloading the tip and 
proplet, while the proplet's effect at high 
loadings benefits the diameter reduction cases. 

Rotational speed reduction is more favor- 
able for the propleted propeller than the 
straight blade propeller. As the rotational 
speed is reduced to decrease the noise, the 
loading increases. Since proplets perform bet- 
ter at the higher loadings, there is less effi- 
c iency loss . 

Figure 10 presents the trade-offs for a bi- 
blade propeller. Again, the loading shift and 
diameter reduction methods have unfavorable 
trends. Blade number increases are for a fixed 
planform, to support the tip device, and so 
increase the viscous drag while decreasing the 
effectiveness of the tip device. Rotational 
speed reduction is more favorable than for 
either the straight blade or propleted propel- 
ler. This is a result of the more effective 
tip device retaining more of the efficiency as 
the loading increases. Using a rotational speed 
reduction, the three way bi-blade propeller has 
about 1 dBA less noise at the same efficiency as 
the six way constant planform straight blade 
propeller. 

The most interesting curve is for the 
effect of sweep on a bi-blade propeller. After 
an initial decrease in efficiency with decreas- 
ing noise, the trend flattens. If the curve is 
extrapolated to the efficiency of a three-way 
constant planform straight blade propeller, the 
three-way bi-blade propeller would have a noise 
level that is 14 dBA lower than for the straight 
blade propeller. 

SUMMARY 

A unique preliminary design tool has been 
developed with the ability to simultaneously 
optimize the performance and noise of a propel- 
ler. By using a vortex lattice aerodynamic 
analysis, a subsonic compact source noise anal- 
ysis, and an optimization scheme, an effective 
and fast preliminary design tool was created. 

This program was then used to investigate the 
noise/performance trade-offs for a straight, a 
propleted, and a bi-blade propeller. 

While the results demonstrate the utility 
of such a design tool, more importantly, the 
relative merits of the different noise reduction 
techniques are identified for optimized geome- 
tries. In short, for straight blades, increas- 
ing the blade number will increase efficiency 
and decrease the noise. Reducing the rotational 
speed will also reduce the noise but will 
decrease the efficiency slightly. For a pro- 
pleted blade, the most favorable method is to 
reduce the rotational speed, but again there is 
an efficiency penalty for noise reduction. 


Finally, for a bi-blade propeller, either a 
rotational speed reduction or sweeping the blade 
will reduce noise effectively. But, with its 
initially greater efficiency, there is a region 
where the bi-blade propeller has a higher effi- 
ciency and a lower noise level than a similar 
straight blade propeller. 
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Figure 2. - Comparison of analytical and experimental 
thrust coefficients for the LT10282+4 propeller. 
Experimental data from reference [91 




ADVANCE RATIO, J = V/nD 


Figure 3. - Comparison of analytical and experimental 
power coefficients for the LT10282+4 propeller. 
Experimental data from reference E9l 
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Figure 4. - Comparison of analytical and experimental 
thrust coefficients for the 1C160 propeller. Experi- 
mental data from reference [10]. 



. 7X10' 1 



Figure 5. - Comparison of analytical and experimental 
power coefficients for the 1C160 propeller. Experi- 
mental data from reference [10] reduced in value by 
half. 
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Figure 6. - Noise source locations. 
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Figure 7. - Analytical and experimental comparison of 
near-field noise for propeller 1C160 in the plane of 
rotation. Experimental data and theory from refer- 
ences tip, 13-15 ], 
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OASPL, dBA 

Figure 8. - Optimized noise/performance trade-offs 
starting from the three way straight blade (S3) 
baseline configuration. Design point; J = 0. 8, 

C p = 0.134. 
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Figure 9. - Optimized noise/perfonr.ance trade-offs 
starting from the three way propleted (P3) baseline 
configuration. Design point; J = 0.8, C p = 0.154. 



SPL, dBA 

Figure 10. - Optimized noise/performance trade-offs 
starting from the three way symmetrically loaded 
bi-blade (SB3) baseline configuration. Design point- 
J = 0. 8, C p = 0. 154. 
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